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Abstract 



X 

5^ , We present analytical next-to-leading order results for the correlator of bary- 

onic currents at the three-loop level with one finite mass quark. We obtain the 
massless and the HQET limits of the correlator as particular cases from the 
general formula, we also give explicit expressions for the moments of the spec- 
tral density. Calculations have been performed with an extensive use of the 
symbolic manipulation programs MATHEMATICA and REDUCE. 
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Baryons form a rich family of particles which has been experimentally studied 
with high accuracy ]l||. A theoretical analysis of these experimental data gives a lot 
of information about the structure of QCD and the numerical values of its parameters. 
The hypothetical limit N^ —>■ oo for the number N^ of colours which is a very pow- 
erful tool for investigating the general properties of gauge interactions was especially 
successful for baryons [@] . The spectrum of baryons is contained in the correlator of 
two baryonic currents and the spectral density associated with it. To leading order 
the correlator is given by a product of Nc fermionic propagators. The diagrams of 
this topology have recently been studied in detail |^, ^, ^ |^, |^ |^. They are rather 
frequently used in phenomenological applications 0. With the advent of new ac- 
celerators and detectors many properties of baryons containing a heavy quark have 
been experimentally measured in recent years [|l|. However, theoretical calculations 
beyond the leading order have not been done for many interesting cases. In this note 
we fill up this gap. 

We report on the results of calculating the a^ corrections to the correlator of 
two baryonic currents with one finite mass quark and two massless quarks. We give 
analytical results and discuss the magnitude of the a^ corrections. The massless and 
HQET limits are obtained as special cases. We also present analytical results for 
the moments of the spectral density associated with the correlator. The extensive 
discussion of the impact of our new result for the correlator on the phenomenology 
of baryons will be given elsewhere. Note that the massless case has been known since 
long ago [T^. The mesonic analogue of our baryonic calculation was completed some 
time ago [|lll] and has subsequently provided a rich source of inspiration for many 
applications in meson physics. 

A generic baryonic current has the form 



j = e'^''^«C4)rvi/, 



which has the quantum numbers of a J^ = 1/2^ baryon for F = 75. \1/ is a finite 
mass quark field with the mass parameter m, u and d are massless quark fields, C is 
the charge conjugation matrix, e"'^'^ is the totally antisymmetric tensor and a, 6, c are 
colour indices for the SU{'i) colour group. Other baryonic currents with any given 
specified quantum numbers are obtained from the current in Eq. (0) by inserting the 
appropriate Dirac matrices. Such additions introduce no principal complication into 
our method of calculation. In the following we take F = 1. The correlator of two 
baryonic currents is expanded as 

ij{Tj{x)3me'^^dx = 7.g^n,(g2) +^n„(g2). (2) 

F = 75 leads to the trivial change 'n.q{q^) — > —Tlq{q'^). The result for the invariant 
function limio^) has already been presented in Ref. |l^. In this note we show results 
for the function ng(g^) and compare it with limiff)- The invariant functions liq^raio^) 
can be represented compactly via the dispersion relation 

where p*(s) = p'''™(s) are the spectral densities. All quantities are understood to 
be appropriately regularized. The spectral density is the real object of interest for 
phenomenological applications and we limit our subsequent discussion therefore to 
the spectral density 

p*{s) = s^ [ptis) (1 + f In (^)) + %f (.)} . (4) 

Here p is the renormalization scale parameter, m is a pole mass of the heavy quark 
(see e.g. Ref. [1^) and a^ = as(/i). The leading order two- loop contribution is 
shown in Fig. |l|(a). Note that this topology coincides with water melon diagrams 
for which a general method of calculation (with arbitrary masses) has recently been 
developed [|, |, 0]. The leading order results read 

pUs) = ^-2z + 2z'- \z^ - 3^2 In z, (5) 



(a) 
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Figure 1: The calculated (a) two-loop and (b-e) three- loop topologies 

pI^(s) = 1 + 9z - 9z^ - z^ + 6z{l + z)\nz (6) 

with z = rn^/s. The next-to-leading order contribution is given by three-loop dia- 
grams with one external momentum. For an arbitrary mass arrangement such dia- 
grams have not yet been calculated analytically. However, if we take the case of one 
massive line, the result within MS-scheme can be obtained analytically and reads 
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where Li„(z) are polylogarithms and C(^) is Riemann's zeta function, C,{2) 



TX' 



76. 



The contributing three- loop diagrams are shown in Figs. |l|(b) to (e). They have 
been evaluated using advanced algebraic methods for multi-loop calculations along 



the lines decribed in Refs. |y, |rT|. This result should be compared to 
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Our method of integration is a completely algebraic one and therefore symbolic ma- 
nipulation programs can be used for performing the long calculations. Two indepen- 
dent calculations of some steps were done using REDUCE and MATHEMATICA. 



REDUCE is rather actively used for high energy calculations (see e.g. Ref. [|T4|)- All 
diagrams have first been reduced to scalar prototypes. The integrals over massless 
loops have been performed (for recurrent integration where possible) and one is left 
with the basic integral 

Via, (3; gV-^) = / (^2 + Jy^ _ j^y, 0) 



which is a generalization of the standard object G{a, /3) of the massless calculation |jT5, 



T6| . The integral V is known analytically and it suffices in order to calculate the 
diagrams in Figs. ||(c) and (d). For the calculation of the diagram shown in Fig. 1(e) 
the basic integral V enters as a subdiagram. This subdiagram then is represented in 
terms of a dispersion integral which makes the whole diagram computable in terms 
of the same V with the argument depending on the loop momentum. The final step 
is a finite range (convolution type) integration over this internal momentum with a 
spectral density of the basic integral V. The reduction to scalar prototypes of the 
diagram shown in Fig. |^(e) leads also to a new irreducible block (i.e. a prototype not 
expressible in terms of V") which is related to a two-loop master (fish) diagram. The 



result for this diagram is taken from Ref. [17 



The results given in Eqs. (0) and (H) represent the full next-to-leading order solu- 
tion. Since the anomalous dimension of the current in Eq. (|l]) is known up to two-loop 
order flSl, the results shown in Eqs. (|3) and (||) complete the ingredients necessary 



for an analysis of the correlator in Eq. (^ within operator product expansion at the 
next-to-leading order level. 

We now turn to the analysis of Eq. (|^ and contrast it with the corresponding 
results for p™(s) as published in Ref. [jl2|. Two limiting cases of general interest are 



the near-threshold and the high energy asymptotics. With our result given in Eq. (0) 
both limits can be taken explicitly. The asymptotic expressions can be also obtained 
in the framework of effective theories which can be viewed as special devices for such 
calculations. 

In the high energy (or, equivalently, small mass) limit z ^ the corrections read 

pUs) = ^^ + hnz-^z-6zlnz + 0{z'), (10) 

p™(s) = 9 + 83z-4:7r^z + 2\nz + 50z\nz + 12z\n^z-24:z({3) + O{z'^). (11) 

Therefore we obtain 

P^is) = pLssless(s) (12) 



mp^is) = mj^{f,)pZssUs) = mj^{p)s' |l + ^ (^2 In (^^^ + |^ | (13) 

where Pmassiess("^) is the result of calculating the correlator in the effective theory of 
massless quarks. For the momentum part p'^{s) we retain the 0{m?) correction. The 
relation between the pole mass m and the MS mass rn-^{p) we have used reads 

™ = m^(,){l + ^(lu(£;)+l)}. (14) 

Note that the massless effective theory cannot reproduce the mass singularities (terms 
like 2; ln(2;) in Eq. (pi])). These singularities can be parametrized by condensates 
of local operators. The first m'^ correction in Eqs. (^ and (11) (or, equivalently. 



the m^ correction to the expression in Eq. (|T^) can be found if the perturba- 
tive value of the heavy quark condensate (0|^\E'|0) taken from the full theory is 
added [|l^. The composite operator (^\E') should be understood within a mass inde- 
pendent renormalization scheme such as the MS-scheme. This value (perturbatively, 
(0|^\1'|0) ~ r/i^ ln(/i^/m^)) cannot be computed within the effective theory of mass- 
less quarks. It provides the proper matching between the perturbative expressions for 
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the correlators of the full (massive) and effective (massless) theories. This matching 
procedure allows one to restore higher order terms of the mass expansion in the full 
theory from the effective massless theory with the mass term treated as a perturba- 
tion pOf. The account for the mass term as a perturbation in a massless theory is 
justified at high energies and greatly simplifies the calculations (see e.g. Ref. |]2l|| ). 



Note that the correction of order m? / s to p'^{s) can actually be found in this manner 
because it depends only on one local operator (^\1/) and, therefore, the calculation 
is technically feasible. In case of the function p'^{s) the situation is different because 
there is no gauge invariant operator of dimension two in the effective massless theory. 
Therefore, the mass singularities of the form m^ log(m^/s) should not appear in the 
expansion for p'^(s) at large energies. The result in Eq. (|12D shows this explicitly. 
Note that the absence of such singularities is one of our checks of the correctness of 
the calculation. 

In the near-threshold limit E ^ with s = {m + E^ one explicitly obtains 



.^(^J:^*^^,,J^-M^n{^ 
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The invariant function p^{s) suffices to determine the complete leading HQET be- 
haviour since one has ^p'' + mp^ ^> ii) + 1)phqet for the leading term. We explictly 
checked this relation. In this region the appropriate device to compute the limit of 
the correlator is HQET (see e.g. Refs. [^2|, ^). Writing 

mp^^m, E) = C{m/p, ««) Vhqet(^, p) (16) 

we obtain the known result for Puqet{E , p) [Q 

,^ , 16E^ (^ as f 182 47r2 ^, a \] ^^.^ 



with the matching coefficient C{m/ij,, as) given by |^| 

C(»/,,„.).l + ^(l,„g)-?). (18) 

The matching procedure allows one to restore the near-threshold limit of the full 
correlator starting from the simpler effective theory near the threshold [^. 

Note that the higher order corrections in E/m to Eq. (^) can easily be obtained 
from the explicit result given in Eq. (|^). Indeed, the next-to-leading order corrections 
in low energy expansion read 
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(19) 
(20) 



To obtain this result starting from HQET is a more difficult task requiring the analysis 
of contributions of higher dimension operators. 



00 

'c£20 r 
2l5 - 



exact solution 
threshold limit 
massless limit 



a 



10 



5 








-5 







2 



4 



12 14 16 
s (in units of m^) 



Figure 2: The ratio Pi" {s) / p^ {s) of the next-to-leading correction and the leading 
order term in dependence of the energy square s 



We now discuss some quantitative features of the correction given in Eq. (|^). 
Of interest is whether the two limiting expressions (the massless limit expression as 



given in Eq. ( p!2[ ) and the HQET limit expression in Eqs. (^) and (|T6D) can be used 
to characterise the full function for all energies. 

For this discussion we compare components of the baryonic spectral function in 
leading and next-to-leading order. In Fig. ^ and ^ we show the ratio Pi{s)/pq{s) for 
# = m and # = q', respectively. In the following we shall always use the specific 
renormalization scale value /x = m if it is not written explicitly. One can see that a 
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Figure 3: The ratio Pi{s)/pq{s) of the next-to-leading correction and the leading order 
term in dependence of the energy square s 



simple interpolation between the two limits can give a rather good approximation for 
the next-to-leading order correction in the complete region of s. 

The informative set of observables are moments of the spectral density 
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where M^ are dimensionless quantities. We find 
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where 



and 
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(n + l)n(n-l)2(n-2)(n-3)' ^ ' 

12 
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S* = A* + ^. (25) 



The coefficients A* are rational numbers, and the closed form expressions for 5* are 
long. Therefore we only show the first values for A^ in the second column of Table |I[ 
They are numerically rather close to the results obtained for A^ in Ref. [|l^ . 
By representing the moments through 



M* M^ 



# 



M#(°) M#(°) 



{l + ^(^*-4)} (26) 



we have all corrections to be normalized to the moment M* of fixed order N. Note 
that the difference S^—S* is scheme-independent. This feature can be used in the high 
precision analysis of heavy quark properties ||^ within NRQCD (see e.g. Ref. P5[] ). 
With Eq. ( |26D one can find the actual (invariant or scheme-independent) magnitude of 
the correction. Indeed, for any given A^ one can find a set of perturbatively commen- 
surate moments M* with n ~ A^ for which the requirement of the chosen precision is 
satisfied. In the third column of Table |I| we therefore present numerical values only 
for the differences of the 61 for subsequent orders. 

Note that the moments represent massive vacuum bubbles, i.e. diagrams with 
massive lines without external momenta. These diagrams have been comprehensively 
analyzed in Refs. p9| , |30| . The analytical results for the ffist few moments at three- 



loop level can be checked independently with existing computer programs for symbolic 
calculations in high energy physics (see e.g. Ref. pT|). 



The presented results have phenomenological applications within the sum rule 
analysis of baryon properties (see e.g. Refs. [0, |33|, ^). As an example we calculate 
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the integral of p^(s) up to the energy cut y/so, 

Mliso) = / ^ p'{s)ds (27) 

which is related to the coupling constant (residue) of a baryon to the current in 
Eq. (|l]). In NLO the integral is represented by 

Mliso) = Mf\so) (l + ^ (in (^^) + A(so))) (28) 

which leads to the renormalization of the LO result for the residue in the form 






Z'n' = ^^^^^ = l + ^ ^- B, +Aiso) +Oial). (29) 



For e.g. So = 2m'^, /i^ = m'^ we find numerically 



^01 = 1 + — A(2m2) = 1 + ^15.4117 . . . (30) 

TT 77 

We see that the NLO correction to the residue in the MS-scheme is rather large. 
For the numerical value of the coupling constant a^ ~ 0.3 which is a typical value 
for baryons containing a c-quark, the NLO correction in the MS-scheme reaches the 
100% level. 

One can see that the corrections to the moments basically reflect the shape of 
the correction to the spectrum. Even the massless approximation is reasonably good 
for relative corrections for the first few moments despite of the unfavorable shape of 
the weight function 1/s". It can be improved by changing the subtraction point /i, 
i.e. by switching from the MS-scheme to some other renormalization scheme, or by 
resumming the integrand |^ which lies beyond the scope of finite order perturbation 



theory though. 

To conclude, we have computed the next-to-leading perturbative corrections to 
the finite mass baryon correlator at three-loop order. Technically, the method allows 
one to obtain analytical results for two-point correlators of composite operators with 
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one finite mass particle that can be compared to HQET results. Corrections in 
E/m near threshold are easily available from our explicit results. From threshold to 
high energies the exact spectral density interpolates nicely between the leading order 
HQET result close to threshold and the asymptotic mass zero result. Going even 
one order higher it is very likely that the full four-loop spectral density can be well 
approximated by the corresponding massless four-loop result which can be calculated 
using existing computational algorithms [^, |37|| . 
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n 




51 - SLi 


4 


9/2 




5 


22/3 


2.833333 


6 


109/12 


1.750000 


7 


5593/540 


1.274074 


8 


6133/540 


1.000000 


9 


460351/37800 


0.821190 


10 


40553/3150 


0.695370 


11 


148574/11025 


0.602132 


12 


2470739/176400 


0.530357 


13 


758614613/52390800 


0.473463 


14 


156200257/10478160 


0.427302 


15 


4583939335/299675376 


0.389128 


16 


117273501721/7491884400 


0.357055 


17 


3113341968041/194788994400 


0.329746 


18 


3172990981751/194788994400 


0.306224 


19 


54887116886639/3311412904800 


0.285760 


20 


11 1547839702373/6622825809600 


0.267802 


21 


7313770708819951/427834547300160 


0.251920 


22 


1483100149208267/85566909460032 


0.237779 


23 


142724395992842749/8128856398703040 


0.225109 



Table 1: Values for the rational part A^ of the first moments 6!^ and their relative 
difference (5^ — S^-i 
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